for i=1:5:300
    for j=1:300
        f=@(x)(i*x+1025)/(x^2+i);
        x0=f(x0);
        fprintf('%i,%.19f\n',i,x0);
    end
end
    